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Abstract 

We study monopole operators at the conformal critical point of the theory in 

2 + 1 spacetime dimensions. Using the state-operator correspondence and a saddle point 
approximation, we compute the scaling dimensions of these operators to next-to-leading 
order in 1/Nb. We find remarkable agreement between our results and numerical studies of 
quantum antiferromagnets on two-dimensional lattices with SU{Nb) global symmetry, using 
the mapping of the monopole operators to valence bond solid order parameters of the lattice 
ant iferromagnet. 
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1 Introduction 


In 2 + 1 dimensions, pure U(l) gauge theory confines (^. One can prevent confinement by 
introducing a sufficiently large number N of massless matter fields, in which case the infrared 
dynamics is believed to be governed by a non-trivial interacting conformal field theory (CFT). 
Such CFTs arise quite frequently in the description of quantum critical points of condensed 
matter systems in two spatial dimensions [2 15 . They also serve as useful toy-models for 
more intricate four-dimensional dynamics, as they can be studied perturbatively in the 1/N 


expansion, where the gauge interactions are suppressed 16 -19 


Our goal in this paper is to study monopole operators in one such CFT, namely the 
theory tuned to criticality. This theory is a nonlinear sigma-model with 
target space, and can be equivalently described as a U(l) gauge theory coupled to Nf, complex 
scalars of unit charge that satisfy a length constraint; see for a textbook treatment. The 
action is 


5 = ^ 
0 


(TX 




1 ) 


( 1 . 1 ) 


where a = 1... iVf,, A is a Lagrange multiplier imposing the length constraint, and g is a 
coupling constant. This theory becomes critical provided that one tunes the coupling to 
0 = 0c for some 0c. 

The interest in monopole operators in this theory is motivated by their interpretation 
as order parameters for the valence bond solid (VBS) order of quantum antiferromag- 
nets 21-23 . The quantum antiferromagnets are defined on bipartite lattices in two spatial 


dimensions, and have a global SU {Nf,) symmetry. Each site of the first (second) sublattice has 
states transforming under the fundamental (anti-fundamental) of SU{Ni,). The sites interact 
via short-range exchange interactions with SU {Nf,) symmetry. There is no explicit reference 
to a gauge field in the lattice Hamiltonian. Nevertheless, when the spin states on each site 
are represented in terms of ‘parton’ degrees of freedom, a U(l) gauge field, emerges in 
the path integral formulation in the 1/Nf, expansion. The partons become the (pa matter 
fields in this gauge theory. As their exchange constants are varied, such antiferromagnets 
can exhibit ground states with two distinct broken symmetries. First, there is the state with 
antiferromagnetic order, in which the SU {Nf,) symmetry is broken by the condensation of 
(pa'- this is the Higgs phase of the U(l) gauge theory. Second, we have the state with VBS 
order in which SU {Nf,) symmetry is preserved but a lattice rotation symmetry is broken. In 
the U(l) gauge theory, this state appears initially as a Coulomb phase; however, the con- 
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finement of the U(l) gauge theory by the proliferation of monopoles leads to the appearance 
of VBS order in the lattice antiferromagnet. This is a consequence of subtle Berry phases 
associated with the monopole tunneling events in the lattice antiferromagnet, which endow 
the monopole operators with non-trivial transformations under lattice symmetry operations, 
identical to those of the VBS order. The quantum phase transition between these two states 
of the lattice antiferromagnet has been argued to be continuous [^[^, and described by the 
theory in (1.1), with monopoles suppressed at the quantum critical point; the critical 
point is therefore ‘deconhned’. 

From the perspective of the theory, this connection to the VBS order of the 

antiferromagnet is powerful because it allows the monopole operators to be expressed as 
simple, local, gauge-invariant operators of the lattice model. Moreover, the couplings of the 
lattice model can be chosen to avoid the ‘sign’ problem of quantum Monte Carlo, and this 
allows efficient studies on large lattices of the CP^'’”^ CFT at the deconhned critical point 


between the Higgs and VBS phases 24-28 . Block et al. 27,28 have obtained the scaling 
dimensions of VBS operators on a number of lattice antiferromagnets, and here we will 
compare their results with the l/N), expansion for the scaling dimensions of the monopole 
operators in the CFT. 

Unlike the lattice antiferromagnet, the monopole operators of the held theory are 


not dehned simply as products of helds that appear in the Lagrangian (1.1). Instead, the 
monopole operators appear as singular boundary conditions that these helds must obey at 


the point where the monopole operator is inserted [^,29,^. One way of dehning monopole 
operators is the following. The model has SU(iV;,)x U(l)top global symmetry. Under 

the action of SU{Nb), the charged scalar helds 0 q, transform in the fundamental represen¬ 
tation. The U(l)top factor in the global symmetry group is a topological symmetry whose 
conserved current is 




( 1 . 2 ) 


The Dirac quantization condition implies that the conserved charge g = f jo satishes 
q E 7^12. It should be noted that our dehnition of q here dihers by a factor of 2 from earlier 


work 23,27,31 


One can dehne the monopole operators as operators that have non-vanishing U(l)top 
charge q. For each q, we will focus on the monopole operator A4q with the lowest scaling 
dimension. The other operators in the same topological charge sector can be thought intu¬ 
itively as products between the monopole operator with lowest scaling dimension and more 
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conventional operators from the g = 0 sector. Throughout this paper, we assume without 
loss of generality that g > 0. The physical quantities that we compute depend only on |g|. 

We aim to determine the scaling dimension of the monopole operator Mg with U(l)top 
charge g to next-to-leading order in 1/W- The most convenient way of performing this com¬ 
putation is to use the state-operator correspondence, under which a local operator inserted at 
the origin of is mapped to a state of the CFT on the conformally-flat background x M. 
The scaling dimension of the monopole operator Mg is therefore mapped to the ground 
state energy on in the sector with magnetic flux J F = dvrg through the 8“^ 29 -31 . (See 


also 23 for an earlier approach to computing scaling dimensions of monopole operators.) As 
is standard in thermodynamics, this ground state energy Ag can be related to the partition 
function on 5^ x M via 


A, = - lim 1 log Z,{I3) = - log Zf x' = 7,, 

/3^oo p ^ 


(1.3) 


where in the middle equality we regularized the 5^ x M partition function by compactifying 
the M direction into a large circle of circumference 13. 

In the case at hand, the 8“^ ground state energy in the presence of 47rg magnetic flux is 
easily computed at leading order in Nb, where the Lagrange multiplier held A and the gauge 
held don’t huctuate and assume a saddle point conhguration that minimizes this energy. It is 
reasonable to assume that the large W saddle point (for both the magnetic hux through 8'^ 
and the value assumed by the Lagrange multiplier held) is rotationally-symmetric. At leading 
order in Nb, the ground state energy on 8“^ comes from performing the Gaussian integral over 


the matter helds 23 . The 1/Nb correction to this result takes into account the Gaussian 
huctuations of the gauge held as well as of those of A. We perform an analysis of these 
huctuations around the rotationally-invariant saddle point; by computing their determinant, 
we extract the 1/Nb correction to the scaling dimension of Mg. Our results are listed in 
Table [T] below. 

In Figure we compare our result for IF 1/2 to numerical studies of the lattice antiferro- 
magnet and hnd a remarkable agreement. This agreement suggests that the next correction 
to J-i /2 in 1/Nb is probably quite small. From comparing the scaling dimensions collected 
in Table to 3, we can also estimate the upper bound on Nb below which the monopole 
operators are expected to be relevant; these bounds are also presented in Table There is 
inherently some uncertainty in these estimates, as they come from extrapolating the large 
Nb expansion to small values of Nb. Nevertheless, our relevance bounds come close to what 


Ref. 27 found from numerics, as can be seen from Table I in 27 
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Q 

A. 

Nb for which Ag < 3 

0 

0 

< oo 

1/2 

0.1245922 Nb + 0.3815 + 0{N-^) 

< 21 

1 

0.3110952 Nb + 0.8745 + OiN-^) 

< 6 

3/2 

0.5440693 Nb + 1.4646 + OiN-^) 

< 2 

2 

0.8157878 Nb + 2.1388 + OiN^^) 

none 

5/2 

1.1214167 Afo + 2.8879 + 

none 


Table 1: Results of the large Nfy expansion of the monopole operator dimensions Ag obtained 
through calculating the ground state energy in the presence of 2q units of magnetic flux 
through S"^. In the last column of the table we listed our estimates for when the monopole 
operators are relevant. 

From Tableit can be noticed that Ag grows with q approximately as In Section]^ 
we give an argument for this behavior, and hnd that 


Ag = q^/^ (0.26408 Nt + c + 0{N-^)) + 0{q^/^) 


(1.4) 


with a constant c ~ 0.23 that can be deduced from the numerical values presented in Table 
In Section we take the first steps towards deriving the value of c analytically. It would be 
very interesting to develop a more complete understanding of the large q behavior. 

Similar calculations were performed in theories with fermionic matter in 32,33 ^ (See 


also 35 where only the fluctuations of the Lagrange multiplier field were calculated with the 


purpose of studying monopole insertions in theories with global U(l) symmetry.) In 36,37 


monopole operators were studied holographically. For studies of monopole operators in 


supersymmetric theories, see, for instance, 38 43 


The rest of this paper is organized as follows. In Section we set up our computation. 
In Sectionwe review the leading order analysis at large Nb. In Section]^ we examine the 
1/Nb corrections around the spherically-symmetric saddle points of the effective action for 
the gauge field and Lagrange multiplier. In Section we analyze the behavior at large q. 
We end with concluding remarks in Section Several technical details of our computation 
are included in the Appendices. 


Whe calculation in the case with scalar matter performed in this paper is technically more challenging 
than that for fermionic matter due to additional UV divergences. One can renormalize these divergences 
using zeta-function techniques, as we do here, but the same result can be derived in other renormalization 


schemes 34 . 
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Figure 1: The scaling dimension of the g = 1/2 monopole operator, J^i/ 2 - The solid line 
is the Nh = oo result (Ref. [^), and the dashed line is the leading 1/Nj, correction com¬ 
puted in the present paper (see Table [^. The quantum Monte Carlo results are for lattice 
antiferromagnets with global SU(A'"b) symmetry on the square (Refs, honeycomb 

(Ref. 


27 


and rectangular (Ref. [27]) lattices. 


2 Setup 


In order to study the large Nj, limit of the CP'^*’ ^ theory, it is convenient to rescale the helds 


such that the action (1.1), appropriately generalized to that on an arbitrary conformally-flat 
space with metric tensor takes the form 


S = - d^x\/g{x) 

S J 


n 




where 72 is the Ricci scalar]^ In this paper we will work on x M, which we parameterize 
by coordinates x = (6*, 0, r). 

The monopole scaling dimension is equal to the ground state energy J^q on x M in the 
presence of a magnetic flux J F = Anq through the Our main task is to determine the 

^We could have absorbed the conformal mass term ^ |baP by shifting iX, but we chose not to. 
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l/iVfo expansion of this ground state energy, which we write as 


T, = Afi (jy + + 0(1/Afft j , 


( 2 . 2 ) 


When q = 0, the corresponding ground state energy Tq is nothing but the scaling dimension 
of the unit operator. We therefore must have = 0- 

It is not hard to see that at large W, the fluctuations of A and of around any saddle 
point configuration are suppressed. Indeed, upon integrating out the scalars in the action 


(2.1), one obtains an effective action for the gauge field and Lagrange multiplier given by 


trlog ( -(V;, - ^ ~ 


iA] — iVfc 

Let’s expand and A around a saddle point by writing 


(2.3) 


iA = yUg + icr, 


(2.4) 


where and a are fluctuations around the saddle point configuration A^ = and iX = jii. 


As can be easily seen from (2.3), the effective action for these fluctuations is proportional 


to W, so their typical size is of order l/i/W and are therefore suppressed at large W- To 
leading order in 1/W, it is therefore correct to set = a = 0, provided that the background 
values Al and are such that the saddle point conditions 


iA] 




iA] 


(T=aa=0 


5X 


= 0 


(2.5) 


(T=au,=0 


are obeyed. One can then develop the 1/W expansion to higher orders by integrating over 


the fluctuations and a using the effective action (2.3). 


In this paper we will focus only on saddles that are rotationally-invariant on and 
translationally-invariant along M. These conditions imply that /i^ is a constant and that, in 
the sector of monopole flux J F = dvrg, the background magnetic field = dA’^ is uniformly 
distributed over S"^-. 


= q sin 9d9 A dcp. 


^In terms of the quantity F introduced in 31 , we have ai = F — . 


( 2 . 6 ) 
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One can choose a gauge where the background gauge potential can be written as 


= g(l — cos 6)d(l). 


(2.7) 


(This expression is well-dehned everywhere away from the South pole aX 9 = vr.) The 


saddle point condition (2.5) is satished provided that the constant /i^ is chosen such that it 


minimizes the value of the effective action evaluated when = a = 0. In other words, the 
equation that determines /i^ is 




( 2 . 8 ) 


This equation depends non-trivially on q, and hence so does /i^. 

In the next section, we calculate the coefficient by simply evaluating at the saddle 
point, while in Section]^ we compute the correction 6J^q from the functional determinant of 
the fluctuations around this saddle point. 


3 Nb = oo theory 

At leading order in N^, one can identify 


= S,nlAl, fijl, 


(3.1) 


evaluated for the value of that solves ( 2 . 8 ), or equivalently at large Nf, 




(3.2) 


and with the coupling g tuned to the critical value 0 = 0 c- In other words. 


1 \ 47r 

= tr log ( -(V^ - iAlf + d + 4 ) “ "^ d 


(3.3) 


Using the fact that the eigenvalues of the gauge-covariant Laplacian on 5^ in the presence 
of magnetic flux dvrg are j{j + 1) — 44,45 , we obtain 









The first term in this expression is divergent and requires regularization. The second term 
is also divergent because, as we explain shortly, the inverse critical coupling 1/Qc diverges 


linearly, and so the second term in (3.4) cancels part of the divergence in the hrst term. To 


be explicit, let us deduce an expression for 1/Qc- The saddle point condition (3.2) at g = 0 
can be written as 


dvr 

0 


j=o 


UJ^ + (j + 1/2)2 + ^2 • 


(3.5) 


The theory is critical when the correlators on x M are those obtained by conformally 
mapping the power-law correlators on For a scalar held 0 q,, this prescription yields a 
conformally coupled scalar on x M, for which po = 0. Hence criticality is achieved when, 
to leading order in N^, we have 0 = 0c and pg = 0, with 


dvr 

0c 


3=0 


a;2 + (j + 1/2)2 


(3.6) 


After substituting (3.6) into (3.4), the resulting expression is still divergent, but can 


be rendered hnite using, for instance, zeta-function regularization as in 35 , or Pauli-Villars 


regularization as in 23,31 . We will not repeat that calculation here. The regularized ground 


state energy coefficient is 


■^r = 2E 


(j + 1/2) [(j +1/2)'^+- ,2]- (j + inf -nn,- f) 


3=q 


2 , ?(1 + V) 

“ + Q 


(3.7) 


This expression can easily be evaluated numerically for any Note that the same procedure 
gives = 0, as required by conformal symmetry. 

As mentioned above, the value of can be obtained from the saddle-point equation 


(2.8), which yields 


E 


3=q 


J + 1/2 


y(j + l/2)2 + /i2_^2 


g = 0. 


(3.8) 


Note that one can obtain this equation directly by differentiating (3.4) with respect to 


without the need of zeta-function regularization. Upon substituting the solution of (3.8) 
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into (3.7), one obtains the values of given in Table These values agree precisely with 


Q 

ha 

'COO 

0 

0 

0 

1/2 

-0.199806 

0.1245922 

1 

-0.397830 

0.3110952 

3/2 

-0.595457 

0.5440693 

2 

-0.792936 

0.8157878 

5/2 

-0.990344 

1.1214167 


Table 2: A few values for the parameters and 
those obtained in 23 by a very different method]^ 


'OO 

9 * 


4 1/A^6 corrections 


In this section we compute the next to leading oder correction to the dimensions of monopole 
operators. The systematics of the calculation are presented in the hrst four subsections, with 


the numerical results presented in Section 4.5 


4.1 Effective action at quadratic level 

To obtain the leading correction bTq to the result of the previous section, one should 
consider the quadratic fluctuations of the gauge held and of the Lagrange multiplier around 


the saddle (2.4). Expanding (2.3) at small a„ and a, one can write the quadratic term in 


the effective action as 


C(2) _ c(2) , c{2) , c(2) 

^(2) _ 

"" 2 , 

d^xd^x \/ g{x) \/ g{x')a{x)D'^{x, x')a{x '), 

5(2) _ ^ 

^aa 2 

d^xd^x'\/g{x) \/ g{x')a^{x)K^'^^ (x, x')a^i{x 

= iV» / 

^ d^xd^x'\/g{x)\/g{x')a{x)H‘^'^'{x, x')a^/{x '), 


(4.1) 


where K'^, and are integration kernels whose expressions will be given shortly. This 
effective action is non-local because it was obtained after integrating out the helds (pa, which 


quantity should be identified with 2p2q in 23 , while our should be identified with —C 2 q 

in 2^. 
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are massless. The kernels appearing in (4.1) can be written in terms of correlators of |0c 
and of the current 




(4.2) 


NtD'’(x,x') = 

NtK'>-'“‘'(x,x') = -\{J'‘(x)r'(x')), + -g‘“‘'S(x-x')(\Mx)\'^).,, (4.3) 

0 0 

= -^{\Mx)fj>^\x'))g, 

0 


where the delta-function contains a factor of 1/ \/g{x) in its definition. The correlators in 
(4.3) are evaluated under the assumption that the gauge field and Lagrange multiplier are 
non-dynamical and fixed at their background values and iX = /i^. The subscript q 

on the angle brackets in the expressions above serves as a reminder of these assumptions. 


Performing the Gaussian integral over and a, we can write the coefficient SJ^q in (2.2) 


as 


= - logdet'M^ 


(4.4) 


where we defined the matrix of kernels 

/ 


M'^{x,x') = 


D'^{x,x') Hl,{x,x') Hl{x,x') ^ 

H‘^'^{x,x') K‘^'^^/{x,x') K'i'^i/{x,x') 

yH'^^{x,x') K^'^i/{x,x') K‘^\'{x,x') j 


(4.5) 


with i = 9,(j) and the primed indices contracting with the index of the held at x'. The prime 


on the determinant in (4.4) means that when computing the functional determinant we 


should ignore the zero eigenvalues that are required to be present due to gauge invariance]^ 
Our goal in the rest of this section is to calculate a regularized version of this determinant, 
thus obtaining 6J^q. 


“For a more detailed treatment of gauge fixing, see 33 . 
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4.2 Eigenvalues of the integration kernels 


One can start evaluating the expressions in (4.3) in terms of the Nf, 
function G'^{x,x') for the complex scalars, which is dehned by 


cx) limit of the Green’s 


{(j)a{x)(j)p{x')) = gdafsC^ix^x') . 


(4.6) 


Performing the required Wick contractions in (4.3), we obtain 


D'^ix^x') = G'^{x,x')G'^*{x,x'), 

Kl^,{x,x') = D^G\x,x')D^,G^*{x,x') - G^*{x,x')D^D^,G\x,x') 

+ D^G’i*{x,x')D^,G^{x,x') - G‘^{x,x')D^D^,G’^*{x,x') (4.7) 

+ 2g^y5{x - x')G’^{x, x ), 

Hl,{x,x') = G\x,x')D^,G'i*{x,x') - G‘^*{x,x')D^,G‘J{x,x') , 


where = d^ — iA^^{x) and D^i = d^i + iA^^{x') denote the gauge-covariant derivatives in 
the presence of the background gauge held. 


In order to calculate the eigenvalues of the matrix of kernels (4.5) required for (4.4), we 


make use of the rotational symmetry and the translational symmetry along M. These 
symmetries imply that the eigenvectors of this matrix are of the form times an ap¬ 

propriate (scalar or vector) spherical harmonic on We will need the usual spherical 
harmonics Yjm{0,(p), as well as the vector harmonics 


^ijm 


vWTi) 


BY 

^ jm 1 


Ak 


(4.8) 


yu = 

Aj{j+ 1)^/9 


where i,k = &, 0, and = 1. We can decompose the Lagrange multiplier huctuation 

a and the gauge held huctuation in terms of these modes. Because, a is a scalar held, 
we only need the usual spherical harmonics for its mode expansion, a,- is also decomposed 
in terms of while Oj is decomposed using the vector harmonics dfj and We will 
refer to the former vector harmonics modes as E modes and to the latter as B modes, as 
they are the S'^ analogs of E and B modes familiar from other contexts: E and B modes 
have vanishing curl and divergence, respectively, and they transform like the E and B held 
under parity. 
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We can expand each of the kernels in Fourier modes as 


^ ^ jm 


— iu){T — T') 


(4.9) 


Ki,(x,x')= f 

^ jm 

du 


—2aj(r—r') 


Ki,(x, x') = I —J2\ Kf''(^)Xi,ira(e, 4>') + •#>') 


jm 


+ m.iJo', 4>') + RT“’{‘^)y..iU0,4>)xi,„(e', <!>') 

du 


^-iui{T-T') 


k;,,(x, x')= I mt'.ao', ‘t’’) + WjJO', 4>') 


^ — iLo(T — T') 


jm 


(4.10) 


H-Ax,x')^ / ‘^Y.^VA)yiAA4')y;jo'A')e 

jm 

du 


— 2Cj(r —r') 




jm 




(4,11) 


and form the matrix of coefficient^ 


M^a;) = 


D’U) 

Hf{ 


^r(^) 


u) 

-Hf*{u) 

Kf^ 


K]-^{u) 

Txq,EB 

j 

(o;) 

-Hf^*{u) 


(c) 


K^,rE 


{-Hf'iA 


(a;) 


j^q,EE 

(^)/ 


(4.12) 


Note that S'^ rotational symmetry implies that these coefficients do not depend on the 
quantum number m. Note that the matrix Mj(a;) is not Hermitian, which can be traced 
to the fact that in the action (2.1) the Lagrange multiplier held A appears multiplied by a 
factor of i. 

The entries of this matrix are related by gauge invariance and CP symmetry. Gauge 


®We are indebted to Nathan Agmon whose work revealed that the original version of this equation 


contained a minus sign error 46 
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invariance of the kernels in position space imply that at separated point^ 


Hi,{x,x') = 0. 


(4.13) 


Plugging in the decompositions (4.10) and (4.11) in these conservation equations, we ob¬ 


tain that the Fourier space kernel M^(ci;) should have the following eigenvectors with zero 
eigenvalue: 

(o| 0 iol -\/Hj + 1)) M’(w) = 0. M’(w)(^ 0 | 0 -iui + 1)) =0. 

(4.14) 

where the pure gauge eigenvector is written in (a | 5, r, E) components, just like Mj(c(;) 
in (4.12). From (4.14) we can express and in terms of Kf 




The second restriction on the entries of Mj(a;) comes from the CP invariance of the 
theory and the monopole background around which we are working. Under CP the modes 
of a and the B modes of transform in the same way, while the r and E modes acquire 


a relative minus sign. Because the effective action (4.1) is invariant under CP, we conclude 


that there is no mixing between a, B and r, E modes. 

These constraints imply that M^(a;) takes a block diagonal form. For j > 0, 


/ 


M^(a;) = 






0 

0 


gr(^) 


0 

0 

Krcoj) 


0 

0 




J ^ > jO+1) j ^ > J 


(4.15) 


This matrix has eigenvalues: 


4 = 




\q T 1) + CJ i^q,TT 

A TP — ; 7 1\. A 

j(j + i) ' 


(4.16) 


^ There are multiple equivalent ways to see that these equations are true. Using the definitions of the 


kernels in terms of correlators (4.3) they are the consequences of the Ward identity V'^U(a:) = 0. Alter¬ 
natively, (4.1) should be zero for a pure gauge configuration = V^a(a:), for arbitrary a{x). Partial 
integration readily gives (4.13). 
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as well as a zero eigenvalue corresponding to a pure gauge mode. 

When j = 0, the harmonics Xjm and yjm are not dehned, so the matrix M^(a;) reduces 
to the 2x2 matrix 


M«(a.) = 




tiri 

Krico)) 


(4.17) 


In addition, the only remaining vector harmonic loo is a constant on and can be gauged 
away. Thus gauge invariance imposes Kq^^{u) = 0, and the only non-vanishing eigenvalue 
is Dl{u). 


We will derive expressions for the entries of the matrices (4.15)-(4.17) shortly. After 


doing so, we can calculate SJ^q from (4.4). It is convenient to subtract SJ^o = 0 from 
The expression we would like to calculate becomes: 


i=o 


det 'M‘1 
det 'MO 


(4.18) 


Using the expression for the eigenvalues from (4.16) and that 
metry, then (4.18) becomes: 


0 by parity sym- 




du 


logilM+i: 

j=i 


Kf 


[OJ] 




(2j + 1) log ■ 


D%u)Kf^{u) + 






(4.19) 


Explicit expressions for the coefficients in (4.15) can be obtained by inverting (4.9 )-(4.11). 


Let us explain how to do so for D'j{u) hrst, and leave the details of how to perform analogous 
computations for the K and H kernels to Appendix For Dj{u) we obtain: 

D%u)27r5{u-u') = [ dV • 


(4.20) 


Since the LHS is independent of m, we can average the RHS over all possible values of 
m. After performing the average, the RHS becomes invariant under performing a combined 

®Note that as the gauge fixing condition is independent of the monopole background, any possible con¬ 
tribution from the Faddeev-Popov ghosts cancels after subtracting the vacuum contribution. 
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rotation in {6, 0) and {O', 0'), so we can take the limit 6' —)■ 0. We can also use that D^{x, x') 
depends only on r — r' to set u' = u and remove the r' integral. The simplihed expression is 


47r 


D%i^) = / <px yjP) lim SD-(x, . (4.21) 


e'-s-o 

r'->0 m=-j 


2j + l 

It is only the m = 0 term that contributes to the sum. Analogous formulas for the K and 


H kernels are given in (A.3). Using explicit formulas for the spherical harmonics, (4.21) can 
be simplihed further to 


D^u) = / d^x \/g{x)Pj{cos0)D''{x, , 


(4.22) 


where by x' = 0 we mean the limit t', O' —)■ 0. Similar expressions (albeit more complicated) 


can be obtained for the other coefficients appearing in (4.15)-(4.17). 


4.3 Kernels at q = 0 


When g = 0, one can obtain closed form formulas for the entries of the matrix (4.15)-(4.17). 


In this case, the Green’s function G^{x,x') can be obtained by conformally mapping the 
one, namely l/( 47 r|x — x'\), and from one can construct position-space expressions for all 
the kernels in (|4.7). The conformal mapping from hat space gives 


G°(x,x') = 


47 r^ 2 (cosh(r — r') — cosy) ’ 
where 7 is the angle between the 2 points on 

cos 7 = cos 0 cos O' + sin 0 sin O' cos(0 — cf)'). 


(4.23) 


(4.24) 


Since G^{x,x') is real, eq. (4.7) implies that H^{x,x') = 0, and consequently H^'^{u)) = 0. 


4.3.1 The D kernel at g = 0 


Plugging (|4.7|) and (|4.23|) into (|4.22|), we obtain 

) pTT 

dr 




Pj{cos 0) 


I sin OdO- . 

0 2(coshr — cos t^) 


(4.25) 
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This expression can be evaluated by first performing the 6 integral, which yields 

D%u) = — / (ire*‘^^Q,(coshr), (4.26) 

J — OO 

where Qj (x) is the Legendre function of the second kind. We can then expand the remaining 
integrand at large r, 




^ n!r(n + j + 3/2) 


(4.27) 


and perform the r integral term by term. The result can be written as 


D%u:) = 


r((j + 1 + iu:)/2) 


4r((j + 2 + ioj)/2) 


(4.28) 


Note that in deriving (4.28) we encountered no divergences in the sums and integrals we 
performed. 


4.3.2 The K kernels at g = 0 

Next, we aim to find an expression for While the expression for that 


follows from (4.10) is UV divergent (as would be its flat space analog), the following difference 
is hnite: 


Ar(^)-A-r(O) = ^ ly £- 1 ] ■ («9) 

It can be checked that 


1 — cos 6 cosh r 
2(cosh r — cos 6*)^ 


= 


-1 


S2 


4 (cosh r — cos 6*) 


(4.30) 


Substituting (4.30) into ( |4.29| ), integrating by parts twice in the sphere directions, and using 
V|27’j(cos0) = —j{j + l)Pj{cos9), one can easily show that 


_ JiJ + 1) nO 


77 --(c) - iLo^-^(O) = 


71 (c). 


(4.31) 
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Similarly, it can be shown that 


kO,bb 


(u) - _ lz,0(0). 




-0,rE^, ,^ _ ,\/i(i + 1) nO 


Ul 


(4.32) 


r^U,r£/ / \ 

A-’ (CJ) = — ZCJ- 




These expressions are consistent with the reqnirements of gange invariance in (4.15). They 


also agree with the flat-space limit expected at large u and j, which was obtained in 14 
Also, it follows from (4.17) that at j = 0 gange invariance reqnires (0) = 0. From (4.15) 


we also know that Taking the cj —)■ 0 limit implies = 0, 

as is hnite. Assnming these relations, we have 


0,EEfr.\ _ 


^ j T 1) pjQ 


Un = 




+ Co, 


t^0.,EE / \ ^ 7 ^ 0 / \ 

{uj) = —DAuj), 


(4.33) 


-j 2 ^ 


K 


0,rEf, , VJU + 1 ) nO 


un = 


-tUJ- 




where the constant Cq remains to be determined; we will see later aronnd (4.57) that Cq = 0. 


4.4 Kernels for general q 

For general q, there is no simple closed form expression for the 0 q, Green’s fnnction. One 
can determine an integral expression for it by hrst expanding the helds (pa in Fonrier modes: 

/ 7 CXD j 

^o.,jm{^)Yqjm{0,(P)e-^‘^\ (4.34) 

j=q m=-j 

where Yq^m are the monopole spherical harmonics introdnced in 
Nfy, the action for the helds (pa becomes 

S^= - ^ ^ ^ + (i + 1/2)^ + hq ~ (f] , (4.35) 

j=q m=-j 


44, 45] . At leading order in 
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from which we can read off 


{uj)^= 27rfl5(ci; - U:')5afi6jji5mm'Gj{uj) , 


(4.36) 


with 


GM ^ 


a;2 + (j + 1/2)2 + /i2 _ q2 ' 


From (4.6), (4.36), (4.34), and (4.37) we can write G‘^{x,x') as 


du 


—iujir—r') 


n 

j=q 

OO 

= 2E 




-Egjlr-r'l 


G,{u) 


3=q 


where in the second line we dehned the polynomial in cosy 




and 0 is a phase factor discnssed in 45 that can be dehned throngh 


e*® cos( 7 / 2 ) = cos(0/2) cos{e'/2) + sin(0/2) sin( 072 ). 


(4.37) 


(4.38) 


(4.39) 


(4.40) 


The angle 7 was dehned in (4.24), and the energy Eqj is 


E,j = /(i +1/2)2 + ^12-52. 


(4.41) 


4.4.1 The D kernel 


Let ns hrst determine Using (4.7), (4.22), and (4.38), we have 


A-jr “ r , _ f,-iEqjl+E^j„)\T\+iujT 

7’(+) = Y E / <i’+\/7ha,,(»)F,.j.(»)F,,j«(»)— yeTb/— 


(4.42) 
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Performing the r integral, we can simplify this expression to 


DUu:) = 




Eqji + Eqj'> 


2j + 1 + i^qj' + Eqj")'^) . 




where 


(4.43) 


Xd(j,/,/')= / smedeEo,,ie)Eq,q,{e)Eq,q,ie). 
Jo 


The 6 integral can be performed analytically, and we have 

■(2j + l)(2/ + l)(2/ + l) 




327r3 


J f f 
0 -q q 


(4.44) 


(4.45) 


We can check that this result equals (4.28) for q = 0 and, for instance, for j = 0 




tanh(7ra;/2) 


2Tr + (.2]'+ ir 


8cj 


Note that the summation in (4.42) is absolutely convergent. 


(4.46) 


4.4.2 The K and El kernels 


Similarly, for the other kernels we can use (4.7), (4.10), (4.11), and (4.38) to obtain 


Hf^u) = Hf{u) = Kf^{u) = Kf^{u) = 0, 
^2 ■ 




u = 


16q7i i 


E 


Eqji + Eqj/^ 


(2j + l)\/j(j + 1) j'J"=q _‘JEqjiEqjii{oj'^ + i^qj' + ^qj"Y) . 


Kri^) = 7T7 
Kf^ioo) = 


27 + 1 ^ 
j',j"=q 


87r2 


— {Eqj! + Eqj»)(u^ + AEqjiEqjii) 
2EqjiEqjii{u‘^ + {Eqji + Eqjii)^) 

iCj{Eqji — Eqjii) 


^D{j,j',j'') + 


‘‘ (2/ + 1) 


AnEqj! ’ 


E 


T^q,EE 

j 


(2j + 7 ', 7 "=, 

87r2 


3 yJ —Q 
oo r 


UJ) = 

Kf^^{u) = 


E 


2EqjiEqji>{uj'^ + {Eqf + Eqj")'^) 

{Eqj' + Eqj'i^ 


2:^(j,/,/), 


(2j + l)j(j + 1) ji jii^g l‘JEqjiEqjii{uj'^ + {Eqji + EqjiiY) _ 


T ( ' ' ''b -L (2j^ + 1) 
e{3,3 ,3 ) ^ AeE„4i ’ 


87r^ 


E 


{Eqji + Eqjii) 


(2j + l)j(j + 1) jijii^q [‘^Eqji Eqjii {u"^ + {Eqji + EqjnY) _ 


2^s(j,j',j") + X] 


'qj 

{‘23' + 1) 

A^Eqji ’ 


(4.47) 
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where 


Xh(j,/,/)= / de sine tan -F'/e)F,j,{e)F^,„{e), 

= I d9 sin9 


>0 

/(/ + l)-/(/ + l)] 
4 


—- 4g2taii2(0/2)F'h(0)F,,v^ 

(4.48) 

A detailed derivation of these formulas is contained in Appendix]^ The quantity 
appearing in (4.48) was given explicitly in ( |4.44 ). Similar explicit expressions for ZnU.f.f) 
and ■,]'') are given in Appendix [Bj 

Note that while the expressions for and above are absolutely 

convergent, those for Kj’^^{u}), and are not and require regularization. 

In order to regularize the latter, it is convenient to hrst compute the quantities 


- K^’^^iO ), (4.49) 


which are free of divergences, and then add back the appropriately regularized values for 
and As argued in the previous subsection, gauge invariance 

implies Kq’'^'^{u) = K^’^^{0) = 0, but does not immediately determine K^’^^{0) denoted by 


Co in (4.33). We can return now to that issue. Let us hrst examine Kq^{uj) 
(|4.47l), we have 


0. From 


kS'^'(O) = 

j '=0 


-( 2 / + 1 ) 

47rCoj' 




j'=0 


AnE^j! 


(4.50) 


(Each sum is divergent individually, but the combined summation is convergent.) Next, we 
can examine Using (4.47) and doing a bit of algebra, we have 


( 0 ) = 


327r 




-2 - 


f + 1/2 f + 1/2 

3 ' + 1/2 j " + 1/2 


i2j' + 1) 

+ 1 / 2 ) 




(4.51) 
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By using the symmetry between the summation in j' and j", this expression can be written 
further as 




Ibvr 


E 


-1 - 


f + 1/2 
f + 1/2 


(2j^ + 1) 

,.^0 4 ^( 2 ' + 1 / 2 ) 




(4.52) 


Summing over j" we can write this expression as 


A'r’“(o) = - - f; +w +1) 

+ ;ep>r(/ + i/2) 


(4.53) 


Both of the individual sums in this expression as well as the combined summation are 
divergent, but gauge invariance dictates that i^/’'^^(0) = 0. This result should be thought 


of as a prescription. It can also be justihed in zeta-function regularization, in which (4.53) 
gives i^r''(0) = C(0,l/2)/(47r) = 0. 


Next, from (4.47), we can also write an expression for 


/T“(o) = ^E 


j'=o 


-fif + l)(2j^ + 1) 

if + 1/2)3 


(2j^ + 1) 

4^(2' + 1/2) 




(4.54) 


Both terms are again divergent, but, using (4.53), we can calculate 


At-(O) - Ay-(O) = i 5: ^ = iz)S(O). 

i'=o ^ ^ 


(4.55) 


Since it'/’'^^(0) = 0, this equation proves that Co = 0 in (4.33). 

We can now provide alternate, appropriately regularized formulas for 
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and that are manifestly convergent. They are: 


Kr^ioj) 


E 


-{Eqj! + Eqji'){u:‘^ + AEqjiEqjii) , ,, ,, 

2j + 1 2Eqq,Eqq„{u:^ + {Eq^, + Eqq,,)^) ^ 


(2j' + 1) 


AnEqj! 


Kq,EE^ 


UJ] = 


E 


8ir“ 


E 


i^qj' + E, 


qj 


(2j + l)j(j + 1) 2EqjiEqjii{u‘^ + {Eqj/ + EqjiiY) 

, ( 2 / + 1 ) 




SnEqji 


+ Cq , 




u = 


E 


8ir= 


E 


(Eqj' + E, 


(4.56) 


qj 


(2j + l)j(j + 1) 2EqjiEqjii{uj‘^ + {Eqj! + EqjiiY) 




+ 


(2j^ + 1) 

SnEqji 


+ C'g ) 


where 


C ^y- 2/ + l y^ 2/ + l 1 

^ SnEg-/ SnEoi/ 47r 

j'=o 


j =q 


^Oj' 


E 


j =q 


f + 1/2 


^(f + 1/2)^ + /i^ - 


-1 -g 


(4.57) 


The expression for K^E'^{u:) was obtained by simply combining the two snmmations in the 


expression in (4.47). The expressions for TfJ’ (cn) and 77/’ (cn) were obtained by snb- 
tracting i7f^^(0) = 0 from the expressions in (4.47). In (4.57) we discover the saddle point 


eqnation for (3.8) (obtained after tnning the conpling g to the critical valne Qc), thns 

Cq = Q. 


One can check from (4.47) and (4.48) that the gange-invariance relations (4.15) are 


obeyed. Snch a check is most simply performed by matching the residnes of the fnnctions of 
uj at their poles at a; = ±i{Eqji + Eqjn) for each f and j", after symmetrization between j' 


and i "—see Appendix A.2 


Now that we have expressions (4.56) we can begin evalnating the order l/Nh corrections 
to the free energy, which is the snbject of the next snbsection. 
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4.5 Numerical Results 


With the regularized formulas for the kernels in hand, we are almost ready to calculate the 


subleading correction bTq to the free energy using (4.19). Let us write this expression as 

1 r rl 

+ (4.58) 

^ ^ _n 


j=0 


where L'j{u) can be read off from (4.19). (See also (C.47)-( C.48| ).) As shown in Appendix [C| 
at large u and j the integrand in this expression behaves as 


LUu) = 


8p2 


n;2 + (j + 1/2)2 


+ ... , 


(4.59) 


thus rendering the integral (4.58) linearly divergent. There are several ways of understanding 


how to regularize this divergence. One way is to use zeta-function regularization to write 

8/^2 


j=0 


a;2 + (j + 1/2)2 


j=0 


1 . . 8/^2 

5E(24 + l)ArT = 4f‘k(0.V2)=0. (4.60) 


Then one can subtract (|4.60|) from (|4.58|), and evaluate 

qW) - 


_ 1 f dU 1 \ 

" 2 y 5? ' 

i=o 




a;2 + (j + l/2)^ 


(4.61) 


instead of (4.58). This expression is no longer linearly divergent. 


Another way of understanding the subtraction in (4.61) is that the critical coupling Qc, 


which was obtained in (3.6) at leading order in receives l/A"b corrections. A similar 


phenomenon was encountered in 14 when computing the thermal free energy at subleading 


order in l/A^"),. Just as in 14 , it can be argued that 


dvr 

0c 


y ^^ 24 + 1 ) 

j=0 


4 ■ 

cn2 + (j + 1/2)2 


+ 0(l/iV,^: 


(4.62) 


The 1/iVfc term in this expression contributes to through the last term in (2.3) precisely 


as the subtraction implemented in (4.61). This expression can be derived rigorously by 
performing a careful renormalization analysis of theory on x M 34 . 


Even after the linear divergence in has been taken care of, this quantity is still 
potentially logarithmically divergent. This logarithmic divergence cancels when using a 
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regularization prescription consistent with conformal symmetry. In practice, we evaluate the 
integral in (4.61) with a symmetric cutoff: 


(j + 1/2)2 + ^2 ^ ^2 _ 


(4.63) 


This can be thought of as preserving rotational invariance on as the high energy modes are 


insensitive to the curvature of the sphere. Given the kernels, (4.56), and the regularization 
described above, we are able to evaluate 5J^q numerically. To obtain good precision, we 


hrst evaluate (4.61) numerically in a region (j + 1/2)^ + < (A')^; then in the region 

(A')2 < (j + 1/2)2 _|_ |;^2 ^ ^2 replace L'j{uj) in (4.61) with the asymptotic expansion 
derived in Appendix C (accurate up to terms of order 0(1/ [{j + ^)2 + a;2]^^^)) and evaluate 
the integral analytically as A —)■ cx). We notice that the result converges very rapidly as we 
increase A'. (In practice. A' = 10 is already sufficiently large.) Our results for 6J^q are given 
in Table See also Table [T] where we collected the results from Tables and together. 


Q 

SRq 

0 

0 

1/2 

0.3815 

1 

0.8745 

3/2 

1.4646 

2 

2.1388 

5/2 

2.8879 


Table 3: The coefficient 6J^q in the large expansion (2.2) of the ground state energy in 
the presence of 2q units of magnetic flux through S'2. 


5 The large q limit 

Whereas the analytic computation of M^(a;) seems to be a hopeless endeavor for hnite q, 
we found that the g —?• oo limit is tractable. The reason for the simplihcation is that this is 
essentially a flat space limit: reintroducing the radius R of we have a strong magnetic held 
5 = g/i?2 on the sphere at large q. 0 quanta move on Landau levels, which are localized on 
1 /v/B = R/y/q distances, hence they don’t feel the effect of the curvature of the sphere. To 
leading order in 1/q, the problem becomes the analysis of the in a constant magnetic 

held in hat space. 

First, we want to calculate We will be more cavalier about divergences than in the 
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rest of the paper, and write, following (3.7), 


(5.1) 


JJ” = + 1) [U + 1/2)" + 

3=Q 
oo 

= + 2n + 1) \2q{n + 1/2) + /i^ + (n + 1/2)^] ^ , 

n=0 

where we introduced n = j — q, and assumed zeta-function regularization as implicit. The 


resulting saddle point equation for /i^ is (3.8): 


E 


2g + 2n + 1 


=0 '^ 2 q'(? 7 - + 1 / 2 ) + /iq + (77. + 1 / 2 )^ 


= 0 . 


(5.2) 


This equation has a solution only provided that we scale correctly with q, namely 


— ‘^QXo + Xi + 0^- 


(5.3) 


Plugging this Ansatz into (5.2) and only keeping the leading terms, we obtain 


o = E 


^ ^/{n + 1/2) + 


Xo 


A / 1 1 

- C 1 2 ’ 2 


(5.4) 


Note that in obtaining this equation we assumed that n ^ q even though we are summing 
over all positive n. This assumption is justihed because the contribution of n > g is higher 
order in 1/g. The constant xo can therefore be obtained as the root of the transcendental 


equation (5.4). Going to one higher order we can determine Xi- 

3C i + Xo) 


Xi — ~Xo + 


which gives the large q expansion for /i^: 


C (§1 ^ + Xo) 


(5.5) 


/i^ -0.39456 g - 0.00456 + O - 

1 \^q 


(5.6) 


As seen from Figure ( |5.6 ) is in excellent agreement with the values of /i^ obtained from 
solving (3.8) at large g|^ 


®In 23 it was noticed that « —2g/5 at large q. 
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Figure 2: A comparison between the values of /i^ found by solving (3.8) (blue points) and 
the large-g analytical approximation (5.6) (solid black line). 


Plugging back into (5.1), we have that at leading order in q 


= 2, ^ = (29^" C (-1, 1 + Xo) 

= 0.26408 9 =’-'^ 


(5.7) 


Quite nicely, this equation can be understood in flat space terms: the Landau levels of a 
massive scalar held are given by 


En = \/2B{n + 1/2) + m2 = — v/2g(nT172y+”2g^ 

rC 


(5.8) 


and have degeneracy J\f = ^ ^ = 2g, giving exactly (5.1) for the free energy (if we set 

R = 1). Note that the scaling of (5.7) follows from hat space dimensional analysis: 


the free energy density is an intensive quantity of mass dimension 3, hence it has to be 
independent of R, and we get R B3I2 

While we leave a full evaluation of Mj(c(;) at large q to future work, as a hrst step we 
now derive large q behavior of Mj(ci;) when a; = 0 and j q. This corresponds in the 
hat space limit to taking the momentum p ^ BR. At cj = 0, we only have to determine 


Hj, and as the rest of the matrix elements vanish—see (4.15). We can 


obtain a closed-form formula for these kernels by taking the explicit expressions for them, 

Q BB 

and expanding for large q. This is quite a tedious task, especially for Kj , where we have 


to expand (B.18) for hxed j, j', j" and large q. The resulting expressions can be summed 


over j', j" analytically using zeta-function regularization. The results are given by the simple 
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expression 
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( ^ 

WjU + i)xo 

0 

o\ 

C (|, i + Xo) 

WjU + i)xo 

2j(j + l)xi 

0 

0 

8 71 

0 

0 

4j(j + l) (Xo + Xi) 

0 


V 0 

0 

0 

0/ 


(5.9) 


j are 0(yg), this 


One curious subtlety is that while the individual terms in Kj 

0{y/q) contribution vanishes upon summation over j', j". The sub leading terms give the 
result in (5.9). This implies that we have to know the saddle point value of to hrst 
subleading order (5.3) and xi makes appearance in the hnal result. 


We obtain the flat space kernel by the replacement j{j + 1) —)• in (5.9). Our expression 
is valid for p BR. It would be an interesting exercise to obtain the full quadratic effective 
action of the model in flat space in a constant magnetic held. 


The matrix in (5.9) has one zero eigenvalue corresponding to the pure gauge mode. The 


nonzero eigenvalues for j = I are (4.16): 


A1 


0.055251 ± 0.023717i 




0.063044 


We can check that (5.10) agrees with the numerical results—see Figure]^ 



(5.10) 


Figure 3: The numerical results for the three eigenvalues, Ag, A+, and Al are plotted against 
the analytic large q value in black. 


lOpor j = 0 the matrix is 2 x 2, and the only nonzero element is Dq( 0), as in (4.17). 






















6 Conclusions 


We have presented here the leading correction to the large W result for the scaling dimension 
of the monopole operator in the CFT. This correction was obtained by computing 

the Gaussian fluctuation determinant of the U(l) gauge held A^, and the Lagrange multiplier 
A, on X M. Computation of higher order terms in the 1 /W expansion appears to be possible 
by the present methods, but will involve considerable effort. 

Our computation now opens the possibility of quantitatively testing the most subtle and 
novel aspects of the theory of deconhned criticality in two-dimensional lattice antifer- 
romagnets. An important feature of this theory is the connection between the monopole 
operator and the VBS operator of the antiferromagnet 2^23 . This connection allows a 
Monte Carlo computation of the monopole scaling dimension by measuring correlators of 
the VBS order in lattice models. We compared our present result with the Monte Carlo 
studies and found very good agreement both in the dimension of the lowest monopole op¬ 
erator and in the number of scalar helds below which monopole operators become relevant; 
see Fig. [^and Table in the Introduction. 

We conclude by noting that our present study is among the most complex theoretical 
computations of critical exponents which have been compared to Monte Carlo simulations. 
Our calculation assumed conformal invariance to realize a framework in which the exponents 
could be determined, and it did not reduce to identifying poles in a Feynman graph expansion 


47 . It would be of great interest to apply the recent progress in bootstrap methods 48,49 


to also determine such exponents. 


Acknowledgments 


We are very grateful to Nathan Agmon, whose joint work with SSP 46 revealed a crucial 


typo in the hrst version of this paper. We thank Ribhu Kaul for valuable discussions and for 
providing the numerical results in Fig. We also thank Max Metlitski for useful discussions. 
The work of ED was supported by the NSF under grant PHY-0756174. The work of MM 
was supported by the Princeton Center for Theoretical Science. SSP was supported in part 
by the US NSF under grant No. PHY-1418069. The work of SS was supported by the NSF 
under grant DMR-1360789 and MURI grant W911NF-14-1-0003 from ARC. Research at 
Perimeter Institute is supported by the Government of Canada through Industry Canada 
and by the Province of Ontario through the Ministry of Research and Innovation. 


29 









A Derivation of integration kernels 


A.l Evaluation of Fourier transforms 


In this Appendix we explain how to obtain the explicit formulas appearing in (4.47). When 
dealing with vectors on 5^ it is convenient to use the frame e^, with z = 0,0 and a = 1,2, 
defined by 


ei = (l,0), e*= 0, 


sin 9 


(A.l) 


One can then convert the coordinate index of a quantity Vi to a frame index by writing 


(A.2) 


The frame index can then be raised and lowered with 5“^ and 6ab, respectively; in other 
words, it makes no difference whether it is upper or lower. 


Our starting point are equations (4.9)-(4.11). In the main text, we explained how (4.9) 


can be inverted to obtain (4.20). The analogous formulas for the other kernels are 


Kr^ioj) 
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Alim V 

m=-j 
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A Mm ^ 

m=-j 

d^xhm Y;je,mHxX)yUo'^ 

va=-j 

J 

d^x lim y W (0, (i))D\x, x')Y.m{0', . 


m=-j 


(A.3) 


Using the spectral decomposition (4.38), the position-space kernels appearing in this expres- 
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sion can be written as 






AEqjiEqji 


^-i^qj'+^qj")\'^-'^'\ 

a;) = Pq,j'i.l)Pq,j"{l) 7? {i^qf ~ ^qj") ~ + ^qj")^i^)] 

^njqjinjqji' 




AnEqj! 


e-i^qj'+^qj")\'r-'r'\ 

Ka'i^^ a:') = J2 - 4 ^ -(^?i' - ^9i") sgn(r - r) 


J J 


X ([D,'e-^‘=^F,j,{y)] [e-‘^'''®F„„(7)]-- [e-“«®F,j.(7)] [£> 76 -“’®F,,7.(7)]*) . 

/c.(^. ^') = E 4^'*“ - ^') + E 


A.En'i^ Ea-i' 

j'J" 


HHx,x') = E 


X ('[-D.e-“«®F,7.(7)] [Ba'e-“’®F,7»(7)]*+ [£>a'e-”’®F,j.(7)] [-0.6““’®-?;,;" (7)]' 
- [e-“«®F,,/(7)] [fl.C<..e-“«®F„,.(7)]' - [£>.£>„.e-“"®F,,7(7)] [e-“«®F,7»(7)]b . 

Q -{ Eqj '+ Eq ^„)\ T - T '\ 


AEn^lEr, 


3 yj 


X [D,e-^^^^EqAl)V-[D.e-^^^^FqAl)] FqAl)V) , 

(A.4) 


where 


Da = da-iAl{x) , Da' =da' +iAl,{x') ■ 


(A.5) 


In taking the x' —)■ 0 limit in (A.3), it is convenient to use the addition formula for the 
spherical harmonics 


E y;j«,mue',r) = -R)j(7) . £ 0 ^( 7 ) = ^^^^(coso), 


dvr 


(A. 6 ) 


m=-j 


where 7 is the relative angle between the points {6, 0) and {6', 0') on S'^. Taking derivatives 
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of this formula and using the dehnition of the vector harmonics in (4.8), we have 


'in E = F„,(e), 

a ;'—>0 
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hm ^ 
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cos (p sm ( 
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0 —Fq-{6)J y— sin0 cos( 


We also have 


lim 

x'^O 


lim Da 

x'—>-0 


lim Da' 




= F,,{e) , 


-igtan ^Fqj{6) j 


x'^O 


lim DaDa' 

x'^0 


= {-F'je) -*,tanfF„(»)) 


cos 0 sm ( 


— sm 0 cos 


0/ ’ (A.8) 




»9 (F,j(e)-sineFC(e)) 
l+cos 9 


-iq(F^j{e)+s\neF'^.{e)) 
l+cos 9 

■^^ + qHw?lFqj{e)^ 


X 


cos 0 sm ( 


— sm 0 cos ( 


Plugging in (A.7) and (A.8) into (A.3) and performing the integrals over oj and 0, we 
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obtain (4.43) and (4.47), with 


where 


/ desmeFo,{e)Fgj,{e)Fg,„{e), 

Jo 

lH(hi\f)= [ d0AaeimJ-Fi^(e)F,y(e)F,y,(0), 

= f dB AnOF^,(e) (F;.,(«)F,j»(») - F„,(«)F;.„(»)) , (A.9) 

Jo 

^EU,j',f)= [ d9sm9AE{9), 

Jo 

[ d9sm9ABi9), 

Jo 


Ae = F” (2F;,f;„ - F„.F''„ - Fd,F„.) - \F„F„,] - ^ tan^ ^F^jF^.F,,,., 


Sin 


2 9d9^-ir^<iri 2 ^°^ 


Ab = 


sin0 


d 


IF,1'F,1''] + Ki (2F;,F;„ - F„,F"„ - Fd.F,j,) 


0 


Vtati“-F|;^F,j.F,j». 


(A.IO) 


To simplify these expressions, we can use the fact that the monopole spherical harmonics 
Yqjjn{9,(j)) are eigenfunctions of the gauge-covariant laplacian on with eigenvalue j{j + 

l)-g2; 

(V^ - - tAl)Yq,U0, 0) = [j{j + 1) - q^] 0) • (A.ll) 

Since Fqj{9) = \J^^Yqj^_q){9,0), we also have 


-^de (sin 980 ) - tan^ ^ + j(j + 1) - 
smfc' 2 


F,,(0)=O, 


(A.12) 


Eq. (A.12) implies that after integration by parts in the formula for above, we have 


MJ,f,f)= / desmeF„j(e)F,j,(0)F,,,„(ff)\j"U" +1) - fu' + 0], (a.is) 


thus obtaining the simplihed form for X^ presented in (4.48). 
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One can simplify Xe as follows. From (A.12), it follows that 


sin^ 


de (sin^F",) = Fh, 


,, , , 1 e 

-J J + 1 + tan - 

sm 6 smO 2 


2q'^ 2 0 / A 1 .x 

+ ^tan -Fgy A. 14 
sm 6 2 


and similarly for Fqjn. After integration by parts in the second and third terms in the hrst 


paranthesis in (A.10) one obtains 


Ae = 2F” + f;, (f'',f;, + f;,f»„) + ^f', {F,,F,,.y 

- f;. [f(f + i)f;.,f„« + j"u" + i)f,,.f;,,] . 


Integrating by parts in the last term and using (A.12), one has 


Ae = 2F''f;,f;,, + f^, (f»,f;„ + f;.,f''„) + ^f;, (f„,f„»)' 

+ ^Oj^qj'^qj” [j'if + 1) + fif + 1)] 




fif + 1 )^+ nf + 1 )^ - ^ [f{f + 1 )+ f{j"+ 1 )] 


COS^ I 


Integrating by parts in the fourth term we get: 


(A.15) 


(A.16) 


+ Foj [KyKy + KyKy) + 


COS^ I 


- K, XrFyid +1) + F;,Fri„j''(j'' +1)] 

- F«y,fF,r U'W + 1 ) - fU" + !))"■ 


Combining the third and fourth terms using (A. 12) gives: 


Ae = 2f»f;,f;„ + f', (f'',f;.„ + f;,f»„) + f;,f;„ 


sin 9 


sin^Fh,)' 


1 


- Fo,F,yF,y, {f{f + + 1 ))^ 

Everything except for the last term is a total derivative, so 


(A.17) 


(A.18) 


Ae = -Fo,F,fF,f, + 1) - f(j" + 1))' 


(A.19) 


The simplihed form for Ze presented in (4.48) immediately follows. The simplihcation of Xe 


is achieved through a similar integration by parts. 
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A.2 Check of gauge invariance 


As discussed around (4.15), gauge invariance implies: 


= y/j{j + l)Kf^{u). 


(A.20) 


We can check that the expression in (4.47) obey these relations by comparing residues at 
u = ±i{Eqji + Eqjii). We have 


Res^=±i(B 


+ {Eqf + Eqjii)‘ 


= T 


Eqji + Eqji, 


From (4.47) and (4.48), one can check that, for instance, when j' ^ j", 


T> (■ {Eqj' 

Res [tujKf {uj}) = I , -,j ) 


2EqjiEqj'. 


2 j + l 

= -Res (\/j(j + , 


(A.21) 


(A.22) 


so the hrst gauge invariance relation is verihed. In checking this relation it was important 
that = {E"^-,, — E‘^-,)ZD{j, j', j”). A similar check shows that the second gauge 


invariance condition in (A.20) is satished. 


B Formulas for Ib and X 


'H 


B.l Results for g = 0, 1/2 


We have only been able to perform the integrals as written in (4.48) for g = 0, 1/2. We had 
to use a different method to obtain results for arbitrary g; the details are given in the next 
subsection. We give the results below. 

For g = 0, we have 




(2j + l)(2j' + l)(2j" + l)(j' + f - j)ij + f - f)ij + f - f) 


327r3 


X 


■ I ■ I ■ ‘ 

3+3 +3 


3 + 3' + f 


'j-l j'-l j"-r 
0 0 0 


(B.l) 
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For q = 1/2, we have 


ll\2 




JU +1) - if - f) 


for j + j' + j” odd 


m + 1) - if + f + 1)^ for J + f + f even 


TBij,f,f') = -XDij,f,f') X 


if-f'?-jij + ^) 


for i + i' + j" odd 


(B.2) 


1 2 


if + f + lf -j U + 1) for j + f + f even 


B.2 Results for arbitrary q 


For arbitrary q we were not able to perform the integrals in (4.48) directly. Rather, we 


employed a method developed in 33 to obtain expressions for the kernels. From these 


expressions we were able to read off explicit formnlas for Xg and Xh- We briefly introduce 
the method below, and give the end results, which are somewhat complicated. 

In this subsection we will use a different basis for monopole vector harmonics than the 
one used in the main text. We denote the harmonics in the new basis by ^q,jmi ^qjm- 

For completeness we provide the details needed for the calculation, but we do not elaborate 


on the technique. For a detailed discussion see 33 


The vector spherical harmonics that we need are tangent vectors to 8 “^ x M. There are 
a few natural choices of basis for this tangent space. The basis we used in the body of the 
text is just the basis induced from the round coordinates on S'^ x M, and is given explicitly 


in (A.l) 


There is also a cartesian basis given by conformally mapping the natural frame basis of 
to X M. We write the standard line element on in spherical coordinates as 

ds ‘^3 = dx'^ = [dr'^ + dO'^ + siif 9d(lf~\ x = ^sin 6* cos 0 sin 0 sin 0 cos 6^ . (B.3) 


The metric on x M is obtained by rescaling the metric (B.3) by e 


(is| 2 xR = d9^ + sifo 9 d(jf + . (B.4) 

We obtain the frame on x M by the conformal transformation of the standard frame 
e ^3 = dx“ on 


e“ = e '^dx° 


(B.5) 
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Explicitly, 


e“ = 


^ cos(0) sin(6*) cos(6*) cos(0) — sin(6*) sin(0) ^ 

sin(6') sin(0) cos(6') sin((;/)) cos(0) sin(6') 
y cos(6*) — sin(6*) 0 y 


The frame satisfies the relation 




ab 


(B.6) 


(B.7) 


There is yet another basis for the tangent space, in which the angular momentum operators, 
J^, and Jz are diagonal. This basis is given by the vielbeins: 


e 


± 



{'^dx + idy) 


= e '^dz , 


(B,8) 


with X, y, z coordinates on 

In this raising and lowering basis the expressions of the vector spherical harmonics for 
j > q are: 


^ j {j-m+l){j-m+2) Y 


Uljmin) = 


(2i+2)(2i+3) 
0-m+l)0-+m+l) T^ 
(i+l)(2i+3) 
(i+m+l)(j+m+2) Y 
(2i+2)(2i+3) 




ijj+l,m+l('h) J 




(- 


2 i 0 +l) \ 

_ m _ Y ■ (n) 


(jr-rrp(j±rn±p-iA 

n 




2j{j+l) ^q,j,m+l[n) j 

Y 2j(2j-l) ^q,j-l,m-l[n) 

/HEybb+yb V (fi\ 

y j{2j-l) 

VY 2j(2j-l) ^q,j-l,m+l[n) J 


(B.9) 


Here, U, V, W all have total angular momentum j, and s runs over {+, 2 :, —}. 

If j = d) we only have two modes U^ j^ifi) and If j = g = 0, only the 

mode is non vanishing. For g > 1, if we have j = g — 1 only the mode Uqjmin) is non¬ 
zero. We don’t encounter this case in this paper. From the raising and lowering basis we 
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can go to the coordinate basis by contracting these expressions with in (B. 8 ). Then 


^q,jm^ ^q,jmy ^q,jm themselves vectors of ordinary monopole harmonics (mnltiplied by 
some coordinate dependent factors coming from the e^). 

A key formnla that we will nse is 




I e-'" 




V^. 


q,jm 


wa 


Q,J^ 


(B.IO) 


where we have to take the scalar prodnct with the coefficient vector. 

After these preparatory steps we can present onr method. We want to calculate the 


matrix elements of the kernel, K.'^Ju), but now in the basis given by (B.9). To do this we 


have to invert (4.10). As explained below (4.20), using rotational symmetry we can average 


over the quantum number m to get a formula that is easier to treat: 


[k;( 


•> m=-j 


(B.ll) 


where X, Z e {U, V, W}. 


}Cj^^,{x,x') is given in terms of Green’s functions in (4.7). We use the hrst line of (4.38) 


to do the Fourier decomposition of G^{x, x') instead of the second line, which we used in the 
rest of the paper. The beneht of treating x and x' separately, rather than introducing the 
relative angle, 7 is an algorithmic method to evaluate the kernels at the expense of having 
long formulas. 


Plugging into (4.7) we get 


r // 


j\m' 

j'\m" 


X 






+ other distribution of derivatives , 

(B.12) 


where we only wrote down explicitly the term coming from the hrst term in (4.7), and we 


grouped spherical harmonics at the same point inside square brackets. For the derivative of 
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the spherical harmonics we can use (B.IO), and plugging this into (B.ll) we get an expression 


involving one scalar and two vector spherical harmonics at the two spacetime points 


We can then use (B.9) to reduce the whole expression to the sum of products involving six 


ordinary spherical harmonics, three at each spacetime point x, As anticipated this is a 
tedious, but algorithmic task. 

We can actually perform the integral over x' with no work. Because we averaged over m 


in (B.ll), the integrand depends only on the relative angle between h and n'. The integral 


with respect to h is therefore independent of h', so we can choose fi' to point in the I direction 
and replace the integral with respect to h' by a factor of dvr. Using 




^q,—m 


2i + l 

Att 


(B.13) 


we get rid of three monopole harmonics. The remaining angular integral over the product 


of three harmonics can be evaluated using some properties of monopole harmonics 45 . 




(B.14) 


and 




= (-!)' 


+ + + (i I' r\ i e 


Att 


a q' q" 


m m m 


(B,15) 


/ j f j" 

where ) is the Wigner 3j-symbol. 

\ m m' m” 

After preforming the frequency integrals over cu', cj", we get the answer as a sum of a 
large number of 3j-symbols, as many of the angular momentum indices on are shifted 


in (B.9). We can use some identities to reduce the number of 3j-symbols that appear in our 


hnal answer, but we are still left with long expressions. The same logic also gives 
and though the latter does not involve any vector harmonics. 


There is an easy translation between the integrals and the kernels (4.47). Because Vo^m 
is identical to the i?-mode of the gauge held, we can read off and Ib from and 

To simplify the expressions somewhat, we introduce the following notation for 


It is more economical to use the frame basis (B.9) to contract vector indices between the vector harmonics. 
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the product of two 3j-symbol that will appear in our formulas 


3 


Af 



3 

3 


' J f 

0 

m 

—m 


0 




\ 0 m 

n 

—n 
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f J f f ^ 

—m J Y 0 —n n J 


(B.16) 


With this notation, we have: 


^hU,3,3‘ 


1 

64 gTT^ 


-(2J + 1)(2/ + 1)(2/ + 1) V(/ -9 + 1) (j' + q) U" -<? + !) ij" + q) X 

< -2(2j + 1)9(2/ + 1)(2/ + 1) X 


+ (2j + 1) (2/+ 1) (2/+ 1) V(j' - q) {]' + g + 1) {]" - q) {]" + 9 + 1) X 


j f f 

0 —q q 
0 q — I 1 — 9 

' ' 3 3" 3' ] 

Q —q q 
(3 q —q 

3 f f 

0 —q q 

0 q + I —q — I 

(B.17) 
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: i2f + 1 ) {-f - j + f + ij "f + f - f i-2q^ + 2 j“ + 2f + 1 )) 


2(22 + l)q^ ( 2 / + 1 ) i2f + 1 ) 

X if + 2f - if + j + 2f + 2,7" - 1) f - if + j + 2f if + 1)) / + f if + 1) i-f - J + r + f)) 
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X (2,7" + l)i-f-j- 2q^ + f + 3 f + 3f + f i-2q^ + 2f + 2f + 3 ) + 2) 


-(2i + 1 ) if + 1 )^ if - g) f if + 1 ) f{-q + j'- 1) (g + f) if - g) (g + f + i) 
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-i2j + l)f iq + f + 1),7" if + 1) i2f + 1) 

X fi-q^-q + f + 3f + 2) if + f - qiq + 1)) i-f -j+ f + f - 2q + f - 2q)f + f + 2) 


(2j + l), 7 '(g' + l)(2/ + l)(g-/'-l)j"^ 

X \/(-g^ + g + f + f) i-q^ + q + f + Sf + 2) {-f -j + if f + f + if + 1) ( 2 g + f + 2 )) 


-2i2j + l)f if + 1 ) (2/ + 1 ) f i-q + f + 1 ) 

X (g + f + 1) i-f - .7 - 2?^ + f - 2 q‘^f + 3f + f i2f + 3) + f ( 2 /' + 3 ) + 2 ) 
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(B.18) 
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One can check that for g = 0, 1/2 we get back the results of the previous subsection. 


C Asymptotic expansions 


The goal of this Appendix is to derive an asymptotic formula at large j and u for the 


integrand in the expression (4.19) for SJ^q. 


C.l Small distance expansion of the Green’s function 


We start by obtaining a better understanding of the scalar Green’s function (4.38). Isolating 


the phase factor appearing in (4.38), we can write 


= e-2*9®Gg(r-r',7), Gq{T - r/y) = —, (C.l) 


OO 


j=Q 




where y and 0 are dehned in (4.24) and (4.40), respectively, and the energy Egj was dehned 


in (4.41). Note that when 9' = (p' = t' = 0, we have Gq{x,0) = Gq{T,6). 
The Green’s function satishes the differential equation 


df + (V^ - zAl){W - - U' + 7 


Gq{T,9) = -:^—^S{T)6{9), (C.2) 

27r sm 9 


where A'^ = q{l — cos9)d(j). It is convenient to change variables to 

. 1 2 ^ • 2 ^ 
t = smh - , s = sm - , 

2 ’ 2 ’ 


(C.3) 


and solve eq. (C.2) at small t and s (where t and s are considered to be of the same order). 


The hrst few terms in this series expansion are 


^ 87ry^Tt ^ dvr 


Aqt + 


s{ — 2Aq + Gq{l + 4/ig))'' 


+ a/s + f 


(g^ + 2/i^)(2s — t) + + t) 


BqE + st 


367r 


+ s" 


2ABq — 4Ag(3 + 4/ig) + Gq(9 + 16g^ + 40/1^ + 16/ig) 

64 


(C.4) 
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where Cg, Ag, and Bg are integration constants. At each integer order in the expansion, 
we have one new integration constant that can be chosen to be the coefficient of t". When 
s = 7 = 0, we have Fgj{0) = {2j + l)/(47r), so 




(C.5) 


The integration constants can be determined from the spectral decomposition (C.l). Indeed, 
for s = 7 = 0, we have 


~ 1 e 


j=<i 


^gj 


(C.6) 


and this expression can be expanded at small r and matched with (C.5). The small r 


expansion is snbtle, however, becanse naively expanding the snmmand at small r resnlts in 
divergent snms. 

In order to only work with absolntely convergent snms, we can hrst consider the qnantity 


^sub _ ^ \ '' 


^-n\T\ 


dTT 

1 

dvr 


n=q-\-l/2 

00 

E 

7i=q-\-l/2 *" 


l + |r|n, 2 2 n . 3 + 3n|r|+n^r^ 


2 n 2 


■(/i -?^) + 






1/2 2\ I 6 2 2\2 I 

len* +■■■ 


(C.7) 


Gg is a power series in fig — q^. The hrst line in this expression was obtained by expanding 


Eqj = ^{j + 1/2)2 _|_ ^2 _ g 2 small fig — q'^ and plngging this expansion into (C.6), thereby 
obtaining an expansion of Gg(0, r). The expression in the second line snbtracts the hrst few 
terms in the small r expansion of the hrst line: at {fig — q^Y order it contains terms np 
to order The snms in (C.7) can be performed analytically, and then expanded at 

small t\ 


^sub _ ^ 

“ 4^ L2v^ 


+ ulVi + 


g2 + 2/i2-3/i; 


1 

+ ■;— 
dvr 


-q + t 


q{l + 2 g 2 - 6/i2) 


9 


+ ... 


+ 


(C.8) 


where in the hrst line we isolated the terms non-analytic in t, and in the second line we 
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Note that this implies that the O {{fig — q^Y) term is absent in the second line. 
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included the terms analytic in t. The non-analytic terms in this expression match exactly 


the non-analytic terms in (C.5), so the difference Gg(0,r) — is analytic in t. In fact, 


using (C.5) and (C.7), it is not hard to see that the difference Gg(0, r) — G^ and its first few 


derivatives w.r.t. r (or t) are expressed as absolutely convergent sums, so one can expand at 
small t (or small r) by expanding the summands. One finds 


G,(0,r) 


oo 


/osub _ 

n=q+l/2 


n 


E( 1 \ 
(^- 2)9 


1 1 +t[2nE(^.^-fil + q^-2n^ 


‘^(nEj -nE^„_i^^-n2(n2-l)-^(3n2-l)(p5-g2)-^(/iJ-g2)2) + 


0-2 5 


(C.9) 


Adding back (C.8), we can obtain an expression for Gg(0,r). It can be checked that this 
expression can be written more succinctly as 


G5(0,r) = - 


1 + 


dTT [2Vi 

1 


9 


dvr 


—^—h Sit + ~ Si)t^ + ... 


(C.IO) 


where the first line contains the terms non-analytic in f, and in the second line Sp is defined 
as the zeta-function regularized sum 


S, = Y.^2j+l){E,jY. 


J=<1 


In particular. 


(C.ll) 




2j + l 
E, 


13 


3=1 

^2 a.^2\ 00 


s. = + E m + 1)B„ - 2U + 1/2)= - + ,=] . 


6 


3=<1 

q{-7 + 10g2 - 18g4 + 30p2(i + 2q‘^) - 90p^) 




240 


E 


J=5 


(C.12) 


(2j + l){E,,r - 2(J + 1/2)^ - 3(/r^ - q^){j + 1/2)^ - - q^ 


q '1 J\J ' / J ^ xr-q 
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and so on. These expressions can easily be evalnated nnmerically. By comparing (C.IO) to 


(C.5), we can extract the constants Cq, Bq, and Aq appearing in (C.IO): 


^ Svr ’ 


A 

^ dvr ’ 


B,= 


S3-S, 

127r 


(C.13) 


Note that the constant Cq is the same as that appearing in (4.57). From (3.7) and (3.8), 
it is easy to see that if one tnnes to criticality we have 


C', = 0, 


'COO 

A = 

dvr 


The constant Bq can be computed numerically for any given q. We have S 1/2 
and Bi ^ -0.01255981. 


(C.14) 


-0.00475009 


With (C.13), we now have a complete expression for the small distance expansion (C.4) 


of the Green’s function. This expansion can be developed to higher orders if needed, and in 


fact the hnal formulas presented below were obtained after keeping one more order in (C.4). 


C.2 UV asymptotic of the scalar kernel 


Let us first examine the scalar kernel D'^{x,x') = \G‘^{x,x')^. From (4.22) and (C.l), we 
have 


D’^Auj) = 2n dr (76^ sin 6* (cos 6')Gq(r, 6*)^ 


(C.15) 


The large j and u behavior of DAui) can be obtained by plugging in the expansion (C.4) in 


(C.15) and evaluating the integrals provided that the terms we’re Fourier transforming are 


non-analytic in s and t. (The analytic terms in s and t are not related in any way to the 
large oj and j behavior of Dj{u).) 

So at large oj and j we have 


DAoj) = 2n j dr dO sin 6 Pj {cos 6) 


P“(5-2t) 


647r2(s-ft) 3271"^ A s + t 

8nBq{3s'^ — 2Ast + 8P) + \32PjjA^ + 4st(—3 + 8/rg) — s^(3 + 20^.^)] 

23 ^ 11 “^ A s + t 


+ 


(C.16) 


In order to evaluate these integrals asymptotically at large u and j, it is convenient to 
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re-expand the term inside the square brackets in terms of simple functions of 


X = iy2(coshr — cos 6*) = 2 a/s + t . 


(C.17) 


For instance, we can write the term inside the square brackets as 


D%u) = F.T. 


T°° R 

" (V|. - 29^) 


{88^ - 248^X12 + 3V^2) 


^OO 


12902407r2 


167r2X2 1927r2 ' ■ '268807r 

-389^ + (29 + 96/iJ)a2v|2 - (23 + 24fil)V%) X^ + ... 


(C.18) 

where we introduced the dehnition of the x M Fourier transform of a function /(0,r) as 


F.T. [/] = 27r / drde sinO Pj {cos 6) f {6, t) . 


(C.19) 


Eqs. (C.18) and (C.16) agree up to the order in s and t to which (C.16) was valid. 


Note that the x M Fourier transform satishes the properties 

F.T. [a;/] = -!.," F.T. I/], 

F.T. [V|,/] =-j(j + l)F.T. [/]. 


(C.20) 


which can be derived upon integration by parts twice in (C.19). These properties, together 


with the explicit Fourier transforms of powers of X given in Appendix C.6 give 

UK.) ^ Z1K-) + T 

^ ^ 2vr [(,+ 1)2+^2]3 ^ [(, + i)2+.;2]^ 

, 8(j + 7)Kl - 3p2) + 3u\j + i)K32/i2 - 23) + 25.^ 


8tt 


[{3 + hr + ‘^4 


+ 0 


[u + iy+<^X\ 

(C.21) 
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Here, D^-{oj) is the quantity defined in (4.28). Asymptotically, at large u and j, 




1 ^ - (j + i) 


8 ^{j + ly + oo^ 64 [ij + 1)2 + n;2] 
ll(j + - 62u\j + 1)2 + lla;4 


5/2 


(C.22) 


1024 [(j + |)2+a;2]®/" 

C.3 The mixed kernel 


+ 0 


[(j + 


7/2 


Next, we examine the mixed kernel H'^[x,x') defined in (4.7). Using G'^{x,x') = e we 

can write 


D'G'^ = e-2*^® [d' + i{A^)' - 2iqd'Q] Gg , 
D'G^* = [d' - i{A^y + 27gd'0] Gg . 


(C.23) 


The definition of then implies 


Hl{x,x') = 2i [Al{A) -2d'^e{x,x')] Gg{T - r', 7 )^ 


(C.24) 


From this expression in can be easily seen that H^{x,x') = 0, which explicitly verifies 
the argument that Hy^{oj) = = 0 because of GP symmetry. For we have 




dvr 


^ 2j + 1 


dr dO dcf) sin 6 lim 

tc'^0 




m=-j 


— ^ [ dr d6 sin 6 e^‘^'^Gl(T, 6)Pj (cos 6) tan- 

vW + i)-/ ^ A 2 


(C.25) 


We should use 


Pf (cos 9) = — sin 9P -(cos 9) 


(C.26) 


and then 



Airiq 

VjU + 1) 


dr d9 sin 9 e*‘^'^G'2(r, 6*)Pj(cos 6*) (1 — cos 9) . 


(C.27) 
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We can integrate by parts and obtain 


+ 1) 


dr dO sin 9 e^‘^'^Pj (cos 9) 


d 


d cos 9 L 


9) (1 — cos 9) 


= (-i)F.T. 


d 


dcos9 


2G'g(r, 9) (1 — cos^) 


(C.28) 


where we nsed the Fonrier transform dehnition in (C.19). We can nse the expansion of the 


Green’s fnnction (C.4), re-expanded in terms of X = -^2(cosh r — cos 9) as in the scalar 
kernel case, to write 




bj] = 


iq 


167r2^j(j -h 1) 


F.T. 


V 52 logX+ — {AdlV%-V%) (XMogX) 




+ ^ - wlv%) x^ + ... 


(C.29) 


This expression agrees with (|C.28|) in a small s and t expansion np to analytic terms in s 

we 


and t. Using the Fonrier transform properties (C.20) and the formnlas in Appendix C.6 


can evalnate the integrals in (C.29) and expand at large u and j. We obtain 




bJ] 


= q 


ij + 


o 


64 [{j + 4)2 +a;2]' 
1 


, 2(j + 4)4_17a;2(j-+l)2+^4 

256 [(j + 4)2+a;2]^ 


(C.30) 


[(j + iy-+ ^ 

C.4 Gauge field kernel 

For the gange field kernel, we have 

^0 = 29 ^Gg(r - r, -f)d^^Gg{T - r, 7) - 2Gg(r - r, j)d^d^>Gg{T - r', 7) 
-4 [A^ix) + 2df,e{x,x')] [Al,{x') - 2d^>e{x,x')] Gq{T-r',-if 
+ 2gf,^A{x - x')G'^{x, x'). 

It is easiest to start with calcnlating the Fonrier modes of dehned by 


(C.31) 


(C.32) 
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Using (C.31), we can extract from 


Kf^"(u) = 2G'iix,x) + 27i dr dO sinO Pj{cos9) 2GgdfGg - 2{drGgf , (C.33) 


where we replaced dr' —>■ —dr when acting on Gg. The first term in (C.33) comes from the 


last line of (C.31). When g = 0, we know from (4.33) that 


-0,rT/ X _ j{j + 1) piO 


^rM= 2 




(C.34) 


where D^{u) was given in (4.28) 


Next, we can calculate the difference K^P^iuj) — K^P^ioj). Integrating by parts the first 


term in (C.33), we have: 


- kY^{oj) = 2 [G'^{x,x) - G°(a;,a;)] 


+ 271 I dr I de sin e P‘^^Pj{cos 6) -AidrGgY - uj^G^ + 4(a,G'o)" + uj^G^ 


(C.35) 


Since from (|C.4[), G^{x,x) — G^{x,x) = Gg, we obtain 


- K+{^) = 2C, - u,^ [B’(a.) - D»(u)] 

— Sir f dr f d9 sin 9 {cos 9) {drGgY — {drGof 


(C.36) 


where we also used (C.15). 


Setting Cg = 0, we can write this expression as 


K(^) - ^?(^)] - F.T. h{drGgY - A{drGof] . (C.37) 
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Expanding at small s and t and keeping only non-analytic terms, we can write 


F.T. 


g2 + 2p2 


= ^F.T. 
Svr^ 


-^IlVs^\ogX + T^{^l-\/l,)X 


irBr, 


480 


(4dlVl, - V‘.) (A'-log A') + ^ (2at - W'tV% + V|,) A 


^OO 

120 


967680 


[39,^ + 5(4pJ - z)dlV% - X^ 


-4(13g" + 20p^)a^V^2 + (64g^ + 92p^)a^V|2 


(C.38) 


(lOg^ + 


{X^\ogX) + ... 


Using (C.20) and the formulas in Appendix C.6 expanded at large uj and j, we obtain 


= Kj’^^{u)+j{j + 1 ) 




+ 


2(j + 


UJ 


_4 [(j + 1/2)2 + ^2]3/2 27r + 1)2 + ^2] 


- 2 g 2 )(^- + 1)2 +( 11 ^ 2 + 


2. . ^ - 27(j + + 4n;^ 


32 [(j + i)2 + c.2] 


7/2 


[(i + i)^ + + ' 


Ar 16(1 - + ir + (88,.J - 87)(J + 1) V - (1 + le^iSo. 


1 ^ 2 , ,2 


,2^, ,4 


Sn 


[ij + iy+u^]' 


(68g2 + 15+2)(j + 1)4 _ 2(478g2 + 111+2)^ 1)2^2 ^ ^§34^2 ^ 5i9;,2)^4 


512 [(4+ 1)2+^2] 


11/2 


+ 0 


[(i + i)^ + + ' 


(C.39) 


Next, we should consider A/’ (a;), which can be computed as 


A/’^^(a;) = 2G"i(+a;) + — 


dvr 


, dr dO dd) sin 9 lim 
2j + 1 J x'^o 




m=-j 


(C.40) 
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where the hrst term comes from the contact term in (C.31). Plugging things in, we obtain 


= 2G’^{x,x) + — 


4:71 


jU + 1) 


dr dO e*' 


dPj {cos 6 ) 

Te 


PUcose)({deG,f -G,dlG, 


Gq ( deGq + 8 q^ esc 6 sin^ -Gq 


(C.41) 


We know from (4.33) that when q = 0, 


f \ ^ P 3 pvO ( \ 

Kj’ M = —-— Dj_^{u), 


(C.42) 


so we can calculate Kf {u) — K^’ {u). Using (C.26) and integrating by parts the hrst 


term on the second line of (C.41), we obtain 


= 2Gq + — 


^ 0 ,BB 


Stt 


Jij +1) 


dr dO e*' 


dPj {cos 9) 

Te 


{deGqf 


d?Pj{cose) ^. 46 -^ dPj{cose) - ^2 


dO^ 


-4q CSC 6* sin -G 


2 ~q 


de 


-{deG, 


(C.43) 


After setting Gq = 0 and integrating by parts once the hrst and third terms under the integral 
sign and twice the second term, we can write 


Kf {u7) = Kf (a;) + 


8vrj(j + 1) 


F.T. 




+ - [4{q^ + + (2/i5 - 3g2)(X^ log X) 

o 'COO 

+ V + ^ [(3 - lQpl)dlV% + {8^^l - 3)V|.] X3 (C.44) 

(17/iJ - 30g2)V|. + 2{8^^l + mq^)dlV% 


40320 


4(16p + 


(X6logX) + ... 
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Using again (C.20) and the formulas in Appendix C.6 expanded at large u and j, we obtain 


'M + 


''COO 


4,/(j + 1)2 + a;2 27r [(j + \f + a;2] 


3(/i^-2g^)(j + |)^ + (9/i^ + 8g>^ 
32 [(j + i)2+a;2]5/2 


+ 6i?o 


(j + 1)^ - 4a;^ 

[ij + lr+u:^]' 


3(28g2 + nl){j + 1)4 _ io(7/i2 + 118g2)(j + l)2a;2 + (808^^ + 459/r2)a; 


512 [(j + l)2+^2] 


9/2 


+ 0 




(C.45) 


C.5 Ultraviolet expansion in the CP^ ^ model 

The ground state energy in the CP'^“^ model can be written as 


j=0 


where we dehned 


Kr^(oj) 


LVu) = log 






when j > 1, and 




LKu,) ^ log 


D§(o,) 


(C.46) 


(C.47) 


(C.48) 


when j = 0. The large j and u behavior of can be easily determined from (C.21), (C.30), 
(C.39), and (C.45) as well as the asymptotic expansion of D^{u) in (C.22). We hnd that for 
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j > 0 we have 




(j + |)2 + a;2+ [(^- + 1)2 +^2] 5/2 

(g^ + 4/i^(8/i^ - 1))(J + 1)^ + 4(-g^ + - 5))a;^ 


12J-- (j + i)2 


■a;" 


144S, 


2[(j + i)2+a;2]5 
3(j + i)^-24(j + i)2a;2 + 8^4 


[(j + ^ 

(25 - 48/i^)(j + 1)^ + 3(64/x2 - 55)(j + + 20a;^ 


(C.49) 


27r 




9/2 


o 


[u+ir+<^^]' 


(From the results of this Appendix, one can construct an asymptotic expansion of L^Aoj) 


that is accurate up to terms that behave as 0(1/ [{j + The derivation of the 


terms not included in (C.49) is straightforward.) The hrst term in (C.49) yields a linear UV 


divergence. We explain how to regularize this UV divergence in Section 4.5 


When j = 0, the analog of (C.49) can be obtained from (C.21) alone: 


700 \ 2 


LS(.) ^ + 32(^r±^ _ 32/^ , 


TTW 


TTbJ'^ 




(C.50) 


C.6 Fourier transforms on S‘^ x 


Here we present some of the Fourier transforms needed in the previous parts of this Appendix. 
Recall that the definition of a Fourier transform of a function of 9 and r on x M is 


F.T. [/] =2ti dr dO sin 6* (cos 6*)/(6*, r). 


(C.51) 


We can calculate explicitly the Fourier transform of 1 7 x 22 ^, where X = ■Y/2(coshr — cos0) 
was defined in (C.17). The calculation proceeds by expanding 1/X2^ at large r and per- 
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forming the integrals term by term: 


F.T. 


X2A 


= 271 drdO sin Oe^'^'^Pj (cos 9) - 


j [2(coshr — cos^)]^ 

2tt\^ [ dr ^^ +i)n(A)j(A — ^/‘^)n ^2n+l -(A+j+2n)\T\+iujT 

(n + 3/2),(27r + l)! 

OO 

2»E 


(C.52) 


(A + j)„(A)j(A l/2)n22n+2 A + J' + 2n 


n=0 


(n + 3/2),(2n + l)! 


(A + j + 2n)2 + oj 


where (a)„ denotes the Pochhammer symbol. The sum in the last expression can be per¬ 
formed analytically in terms of hypergeometric functions, but we will not find it helpful to 
do so here. 

When A = ^ — m with m > 0 an integer, the inhnite sum in (C.52) becomes hnite 
because only the hrst m -|- 1 terms contribute. We have 


F.T. 


1 

X 


dvr 


(j + 5)2 + ’ 


F.T. [X] = 
F.T. [X3] = 


Stt 


Qbvr 


(C.53) 


[u +1)"'+[(j + !)'’+[(j 

and so on. The general formula is 

(-l)"*(2m)!(47r) 


F.T. [X-^+2™] = 


nr.„ [(j +1 - m+ 2k)^+ (.,2] 


(C.54) 


As can be deduced from Section 4.3.1, we have 


F.T. 


X2 


= IQtt^D^Aoj 


(C.55) 


By taking derivatives of (C.52) with respect to A and evaluating the resulting expression at 


54 





















A = —2m, where m > 0 is an integer, we can also calculate 


F.T. [logX] 
F.T. [X^logX] 
F.T. [XMogX] 
F.T. [X®logX] 
F.T. [X®logX] 




_ X 6 _ 0 

[a;2 + (j - 1)2] [e^2 + + 1)2] j'+2l^J , 

_ 16^^ X 120 _^0 . ) 

[a;2 + (^-_2)2][a;2+_^-2][^2 + (_^- + 2)2] ^'+3^ 

levr^ X 5040 / x 

—5--h'T 4 a; , 

nLo K + (J - 3 + 2A.P] '• 

IOtt^ X 362880 

nLo + (i - 4 + 2A;)2] 


DUi^) , 


(C.56) 


and so on. 
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